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ABSTRACT 

AU Microscopii is a 12 Myr old M dwarf that harbors an optically thin, edge-on disk of dust. 
The scattered light surface brightness falls with projected distance b from the star as b~ a ; within 
b = 43 AU, a w 1-2, while outside 43 AU, a w 4-5. We devise a theory to explain this profile. At a 
stellocentric distance r = tbr = 43 AU, we posit a ring of parent bodies on circular orbits: the "birth 
ring," wherein micron-sized grains are born from the collisional attrition of parent bodies. The "inner 
disk" at r < tbr contains grains that migrate inward by corpuscular and Poynting-Robertson (CPR) 
drag. The "outer disk" at r > tbr comprises grains just large enough to remain bound to the star, 
on orbits rendered highly eccentric by stellar wind and radiation pressure. How the vertical optical 
depth tj_ scales with r depends on the fraction of grains that migrate inward by CPR drag without 
suffering a collision. If this fraction is large, the inner disk and birth ring share the same optical 
depth, and tj_ oc r~ 5 / 2 in the outer disk. By contrast, under collision-dominated conditions, the inner 
disk is empty, and t± oc r -3 / 2 outside. These scaling relations, which we derive analytically and 
confirm numerically, are robust against uncertainties in the grain size distribution. By simultaneously 
modeling the surface brightness and thermal spectrum, we break model degeneracies to establish that 
the AU Mic system is collision-dominated, and that its narrow birth ring contains a lunar mass of 
decimeter-sized bodies. The inner disk is devoid of micron-sized grains; the surface brightness at 
b < 43 AU arises from light forward scattered by the birth ring. Inside b = 43 AU, the disk's V — H 
color should not vary with 6; outside, the disk must become bluer as ever smaller grains are probed. 

Subject headings: accretion, accretion disks — celestial mechanics — circumstellar matter — planetary 
systems: formation — stars: individual (AU Mic) — stars: mass loss 



1. INTRODUCTION 

"Debris disks" surrounding young stars are composed 
of optically thin dust (see the reviews by Artymowicz 
2000; Lagrange, Backman, & Artymowicz 2000; Zuck- 
erman 2001). Most debris disks are inferred to exist 
from measurements of infrared excesses (e.g., Aumann et 
al. 1984; Habing et al. 2001; Zuckerman & Song 2004). 
A few disks are close enough to resolve in images, either 
in scattered starlight (e.g., Smith & Terrile 1984; Schnei- 
der et al. 1999; Kalas, Graham, & Clampin 2005) or in 
thermal emission (e.g., Telesco et al. 2000; Greaves et 
al. 2004). 

Is the observed dust primordial — the remains of an op- 
tically thick, gas-rich disk from a previous Herbig Ae or 
T Tauri phase? Or is it maintained in equilibrium — 
continuously removed by processes such as Poynting- 
Robertson drag and replenished by the comminution of 
larger, colliding parent bodies? A third possibility is that 
the observed dust represents the transient aftermath of 
recent cataclysmic events. Dust might be freshly gen- 
erated, uncquilibrated debris from the catastrophic de- 
struction of large planetesimals (Su et al. 2005; Song et 
al. 2005). 

The debris disk encircling the young M dwarf AU Mi- 
croscopii is a promising place to investigate these ques- 
tions. It is well resolved in scattered light from optical 
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to near-infrared wavelengths (Kalas, Liu, & Matthews 
2004; Krist et al. 2005, hereafter K05; Liu 2004; Metchev 
et al. 2005). Of central relevance to our study is the 
disk's surface brightness profile. Within a projected 
distance b from the star of 43 AU, the surface bright- 
ness SB falls approximately as b~ l s (K05). We refer 
to this region as the "inner disk." Outside 43 AU, in 
the "outer disk," the slope of the profile changes dra- 
matically: SB oc b~ 4 - 7 (K05). This break is observed 
independently by other researchers (Liu 2004; Metchev 
et al. 2005). The shape of the profile is all the more 
significant because it r esembles that of the debris dis k 
surrounding /3 Pictoris l|Kalas fc ,Tewittlll995t lLivi r2004) , 
and of the recently disco vered disk encompassing HD 
139664 ijKalas et al J 12006]) . AU Mic's disk is also de- 
tcctcd in unresolved thermal emission (|Liu et alJl200l 
iChen et alJl2005h . The disk's infrared spectrum peaks 
at a wavelength of ^100 /im and exhibits no excess at 
12 (tm; this behavior sugge sts that the disk contains an 
inner hole (jLiu et alJl20"ol . 

Here we offer a theory that explains these observations 
quantitatively. The reason why the surface brightness 
profile breaks at 43 AU is that a narrow ring of parent 
bodies, analogous to the solar system's Kuiper Belt, or- 
bits the star at a stellocentric radius r = tbr = 43 AU. 
The subscript "BR" refers to our term for the belt of 
parent bodies, the "birth ring," wherein micron-sized 
dust grains are born through collisions of larger plan- 
etesimals. Grain creation is balanced in steady state by 
destructive grain-grain collisions and removal by corpus- 
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cular and Poynting-Robertson (CPR) drag. 3 Corpuscu- 
lar drag exerted by the young M dwarf's wind is proba- 
bly at least a few times more effective at removing dust 
than Poynting-Robertson drag in the AU Mic system, a 
possibility first pointed out by Plavchan, Jura, & Lip- 
scy (2005, hereafter PJL05). The outer disk comprises 
grains that are only tenuously bound, moving on orbits 
rendered highly eccentric by stellar wind and radiation 
forces (e.g., Lecavelier des Etangs et al. 1996; Augereau 
et al. 2001). These barely bound grains dominate scatter- 
ing of starlight in the outer disk. By contrast, unbound 
grains escape the system too quickly for their steady- 
state population to contribute appreciably to the surface 
brightness. The inner disk is populated by grains that 
migrate inward by CPR drag quickly enough to evade 
collisional destruction. In CPR-dominated (what we re- 
fer to as "type A") disks, a large fraction of grains meets 
this condition, unlike in collision-dominated ("ty pe B") 
disks. Simi lar classificat i ons w ere put forward bv lWvattl 
(2005) and lMever et al.1 1|2006|) in their considerations of 
disks composed of single-sized grains. We calculate si- 
multaneously the steady-state spatial and size distribu- 
tions of dust particles, and derive how the outer disk's 
optical depth scales with radius for type A and type B 
disks. Our analysis accounts for destructive grain-grain 
collisions and the detailed dynamics of CPR drag, which 
reduces not only the orbital semi-major axes of grains but 
also their orbital eccentricities (Wyatt & Whipple 1950). 
The reduction of orbital eccentricity is not often consid- 
ered but is a key component of our theory. Whether 
type A or type B conditions apply to AU Mic's disk is 
determined in part by the strength of the stellar wind, 
which according to previous works is uncertain by orders 
of magnitude. In this work, we place a novel constraint 
on the stellar mass loss rate and decide the appropri- 
ate disk type by comparing our theoretical models to the 
observations. 

In |21 we lay down basic parameters of the AU Mic sys- 
tem: stellar properties, disk optical depths, timescales for 
grain-grain collisions, and how the star's wind and radi- 
ation alter orbits of dust grains. In J|J we employ order- 
of- magnitude physics and analytic scalings to understand 
how the interplay between collisions, blow-out, and drag 
shapes the observed surface brightness profiles of the in- 
ner and outer disks. There we derive the steady-state 
grain size distribution as a function of position, includ- 
ing the maximum sizes and total mass of parent bodies. 
In U| we verify and extend our analytic results with a 
Monte Carlo simulation of the disk's surface brightness, 
color, and spectral energy distribution (SED). Models 
are compared directly with observations. Finally, in jj^l 
we summarize our theory, place it in context with our 
understanding of how planets form, and point out direc- 
tions for future research. 

2. PRELIMINARIES 

We establish orders of magnitude characterizing the 
AU Mic system. Stellar properties — including the stel- 
lar mass-loss rate that figures prominently throughout 
our analysis — are reviewed in W2.l\ Collision times be- 

3 Our analysis ignores any gas that might still be orbiting the 
star. Only upper limits are observ ed for the column of gas towards 
the star: N H2 < 1.7 X 10 19 cm -2 IRoberge et al.ll2005l) . 



tween grains and relative collision velocities are esti- 
mated in H2.2I Grain dynamics relevant to our theory 
include blow-out by stellar wind and radiation pressure, 
as discussed in H2.3I and orbital decay by corpuscular 
and Poynting-Robertson drag, as treated in H2.4I 

2.1. Stellar Properties 

AU Microscopii is a spectral type dMle star located 
a distance d = 9.9 pc from Earth. It has mass M* — 
O.5M , radius i?* = 0.93i? Q , effective temperature T* = 
3500 K, luminosity L* = 0.1L Q , V magnitude 8.8, and 
H magnitude 4.8 (Kalas et al. 2004; Metchev et al. 2005; 
and references therein). The star's age is estimated to be 
^agc = 12^4 Myr. AU Mic's X-ray luminosity is a prodi- 
gious L x = 5.5 x 10 29 erg/s = 3 x 10~ 3 L* ijHiinsch et alJ 
1999) . Th e star flares at bot h X-ray and ultraviolet wave- 
lengths llMagee et alJl20 031. The stellar rotation period 
is 4.87 davs (|Torres fc Ferraz Mellolll973fl . 

How strong is AU Mic's stellar wind? The wind ve- 
locity v W ind is likely of order the stellar escape velocity 
v csc = ^/2GAU/R* w 450km/s. Plavchan et al. (2005) 
discuss what is known about mass loss rates from M 
dwarfs, citing values ranging from 10 to as high as 10 4 
times the solar mass-loss rate of M Q = 2x 10 _14 M Q yr _1 . 
While the star's youth, flaring activity, and fast rotation 
suggest that a powerful wind emanates from AU Mic, 
the star's X-ray emission indicates otherwise. Wood et 
al. (2002) and Wood et al. (2005) study the relationship 
between X-ray luminosity and stellar mass loss rate by 
measuring M* from a handful of M, K, and G dwarfs 
having ages > 500 Myr. They establish that the mass 
flux at the stellar surface Fm = M*/ (47rF 2 ) increases 
with X-ray surface flux Fx = £x/(47rF 2 ) for Fx < 
8 x 10 5 ergcm~ 2 s _1 = Fx.crit- The mass flux Fm sat- 
urates at a maximum of 10 2 Fm,©. For Fx > Fx lC rit, Fm 
drops to < IOFm,©, perhaps because the strong magnetic 
fields of such extraordinarily X-ray-active stars inhibit 
stellar winds (see also Schrijver & Title 2001; Strassmeier 
2002). AU Mic's X-ray flux F x ~ lx 10 7 erg cm" 2 s" 1 ex- 
ceeds Fx, C rit , implying at face value a relatively low mass- 
loss rate. Nevertheless, it is unclear whether the mea- 
surements of Wood et al. (2005) apply to this highly vari- 
able, pre-main-sequence star. To accommodate the un- 
certainty in AU Mic's mass-loss rate, we consider in our 
analysis a wide range of values, M* £ (1, 10, 10 2 , 1O 3 )M0. 
We ultimately find in 21 that detailed comparison be- 
tween theoretical models of the disk and observations 
can, in fact, constrain M*. 

2.2. Collision Times 

Consider for the moment a disk of single-sized particles 
on low-eccentricity orbits. Where the vertical, geomet- 
ric optical depth equals tj_, the mean free time between 
collisions is 



where is the local orbital angular frequency. This 
expression may be derived by recognizing that every 
~l/fi orbital period, a typical particle traverses an opti- 
cal depth of over the course of its vertical epicycle. 
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For t_l < 1, the particle collides with probability t_l; for 
t_l > 1, it undergoes t± collisions per orbit. 

Detailed fits to observations in 21 reveal that the total 
vertical, geometric optical depth in the birth ring at r = 
r B R = 43 AU equals tx,br = T_L (^br) = 4 x 10~ 3 . We 
define a fiducial collision time 



^c,BR 



1 



■ 2 x 10 4 



4 x 10- 

TL.BR 



yr. 



(2) 



The collision lifetime i co i(s) is the time a grain of radius 
s survives before it is destroyed by colliding with another 
grain. Unlike the case for t c , in calculating t co \ we do not 
assume that particles have a singl e size. For a given colli- 
sional specific energy Q* (ergs g~ 1 lGreenberg et alll978t 
iFuiiwara et alJll989j) . targets of size s suffer catastrophic 
dispersal by smaller projectiles of minimum size 
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2Q* 



1/3 



^0.6 



1/3 



lOOms - 



,10 7 ergg-V V Wrd J (3) 
which is comparable to s. By catastrophic dispersal we 
mean that the mass of the largest postcollision fragment 
is no greater than half the mass of the original target and 
that collision fragments disperse without gravitational 
reassembly. We have normalized the relative speed v re \ 
to the vertical velocity dispersion of visible grains at r = 
r BR = 43AU, 



'BR 



100ms~ 



(4) 



where the full vertical disk height /i BR w 2.75 AU 
(K05). We have normalized Q* to a value appropriate for 
centimeter-sized silic ate targets ijGreenberg et al.lll978t 
IFuiiwara et alJ Il989|) . Ice targets of similar size have 
specific energies that a re 2 orders of magnitude smaller 
ijGreenberg et alll978l) . On the other ha nd, it is possible 
that Q* increases with decreasin g size llFuiiwara et al.l 
119891 IHousen fe Holsappld I1990]) . perhaps as fast as 
Q* cx s~ 5 . If so, grains having sizes s ~ 1 pm would be 
~100 times stronger than their centimeter-sized coun- 
terparts, thereby possibly cancelling the reduction in 
strength due to an icy composition. To keep the telling 
of our story as simple as possible, we adhere to a nomi- 
nal, size-independent value of Q* = 10 7 ergg _1 . The es- 
sential point is that collisions between comparably sized 
grains in the AU Mic disk are destructive. 4 

2.3. Blow-out by Stellar Wind and Radiation Pressure 

Grains of certain sizes cannot occupy orbits bound to 
the star because of stellar wind and radiation (SWR) 
forces. The ratio of repulsive to gravitational forces felt 
by a grain equals 



0-- 



-^rad ~t~ -*wind 
- 1 grav 

3 L^Pswr 

16-7T GM*CpS ' 



(5) 
(6) 



4 Commercial sandblasting machines accelerate abrasive parti- 
cles up to speeds of 100 m s _1 . 



where the dimensionless factor 



-PSWR = Qrad + Qv 



M*V wind C 



(7) 



measures the extent to which the ram pressure exerted by 
the (assumed radial) wind dominates radiation pressure. 5 
Here s and p ~ 2 g cm -3 are, respectively, the radius and 
internal density of a particle, G is the gravitational con- 
stant, c is the speed of light, Q ra d < 2 is the cross section 
that the grain presents to radiation pressure in units o f 
the geometric cross section l)Burns. Lamv. fc Soterl l9791. 
and Q w ind ~ 1 is the analogous dimensionless cross sec- 
tion the grain presents to wind pressure. 6 For Qrad ~ 2 
(appropriate for the s > A*/27r geometric optics limit 
where A* « 1 /xm is the wavelength at which the bulk 
of the stellar luminosity is emitted), wind pressure is 
negligible compared to radiation pressure except for the 
largest value of M* considered. Table Q lists possible 
values for Pswr- 

Grains are continually created from colliding parent 
bodies. Throughout this paper, we assume that parent 
bodies move on nearly circular (/3 <C 1) orbits, and that 
the velocities with which grains are ejected from par- 
ent bodies are small compared to parent body orbital 
velocities. These assumptions imply that grains having 
> 1/2 are "blown out" by SWR pressure. For constant 
-Pswr with s, the condition j3 > 1/2 is equivalent to 



3 £*-Ps 



*-fSWR 



8?r GM*cp 



0.2 



p< 



SWR 



/Ltm . 



(8) 



Grains for which s < Sbiow are unbound and move on hy- 
perbolic escape trajectories. Grains for which s = s^ow 
move on parabolic escape trajectories. A grain for which 
s — s biow "C Sbiow moves initially on a highly elliptical 
orbit whose periastron distance r per i o equals the orbital 
radius of the parent body. It is these "barely bound" 
grains that contribute significantly to the scattered light 
observed in the outer disk. The initial eccentricity eo 
and semi-major axis oq of a barely bound grain upon its 
birth are uniquely related to grain size s via the force 
ratio (3: 



eo 



a 



1-/3' 

^pcri.O 

1 - eo ' 



a _ 1 f fblow\ 
' 2 V 8 ) 



(9a) 
(9b) 
(9c) 



where the last relation assumes that -Pswr is independent 
of s. 

A consequence of equations Q and © is that grains 
on highly eccentric orbits have prolonged collision times. 
Provided that the optical depth traversed by a grain is 

5 The assumption of a purely radial wind is valid insofar as 
the azimuthal velocity of the wind w w i n d,c/> is l e ss than Qr. By 
modelling t he stellar magnetic field as that o f a rotating (split) 
monopole (Weber & Davis 11961 IParkerlll96l) . we estimate that 
«wind,0/£fr £ 1/30. 

6 Our Q r ad equals Q pr from lBurns. Lamv. Soterl I197flfl , It 
should not be confused with Qscat, the usual scattering efficiency. 
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concentrated near periastron at r 



i, the optical 



depth Ti in should be evaluated at r = r pcr i- How- 
ever, £1 should not necessarily be evalu ated for a semi- 
major axis equal to r per i- Instead, from (|9b|l . 



n(e) 



GAU 



1/2 



(1 



13/2 



(10) 



A related useful quantity is the fraction of time a grain on 
a fixed orbit spends at radii between r pcr i and r per i + Ar: 



/(e) 



Ar 



1/2 



(1 



13/2 



■„.,■■/ (1 + e) 1 / 2 ' 
valid in the limit of Ar <C rpcri and large e. 



(11) 



2.4. Corpuscular and Poynting-Robertson Drag 

Under the drag due to corp uscular and Poynting- 
Robe rtson (CPR) forces (see, e.g. . IB-urns. Lamv. fc Soterl 
1979), dust grain orbits having periastron distances r per i 
and arbitrary eccentricities e collapse to a point in a time 



CPR 



3L*Pcpr 
where the dimensionless factor 

-PCPR = Qrad + Q 



£( e ) r pcri S , 



wind " 



(12) 



(13) 



quantifies the relative importance of wind-driven versus 
radiation-driven drag. For the values of M* we consider, 
the stellar wind is at least as important as the stellar 
radiation: Pcpr > 2; see Table ^ The dimensionless 
factor 



E(e) = 



(1 + e) 2 

e 8/5 



c 3 / 5 dx 



(1 - x 2 ) 3 / 2 



(14) 



pverns the decay of orbital eccentricity 
.■Wvatt fc Whipple! 11950ft. A related useful result 
from lWvatt fc Whipple! ljl950|) is that 



de 
lit. 



15L,P C PR 



1 



8/5 



(1 - e 2 ) 3 / 2 



327rc 2 J or 2 



pcri,0 



s (1 +e ) 2 



=3/5 



(15) 



The initial effect of CPR drag on a highly eccentric 
orbit is to reduce the apastron distance while keeping 
the periastr on distance and eccentric ity nearly fixed (see 
Figure 1 of I Wvatt fc Whipple! ll95C|) . For highly eccen- 
tric orbits, the time spent during this apastron reduction 
phase is much longer than the usual CPR-timescales that 
are cited for e <C 1. As e approaches 1, E diverges as 



E(e « 1) oc (1 - e) 



-1/2 



(16) 



Once CPR drag reduces the eccentricity to values less 
than a few percent, the entire orbit collapses to a point 
in a time given by tcPR with E w 1 . 

Table ^ provides sample values of tcpn(e = eo) for 
three choices of s/s D iow = {1.1,1.5,15}, while Figure^ 
depicts how tcPR( e = e o) 

varies with. s, all for rp Cr i — 

r BR = 43AU. 
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Fig. 1. — CPR drag time as a function of grain size s/sbiow: 
for M,/Mq = {1, 10, 10 2 , 10 3 }. To highlight the behavior as 
s — > s blowi the horizontal axis is scaled as log(s/sbiow — !)• The 
fiducial collision time t co \ = tc .BR(^/ ^blow) 1//2 (1 — eo) _3//2 is also 
indicated as a dashed line (see i|3.1.3l . The size Sbreak corresponds 
approximately to where this collision time and the CPR drag time 
are equal. Timescales for removal by CPR drag and collisions can 
be much shorter than the age of the system, t age . 



3. THEORY 

We assemble the ingredients laid out in Shinto a theory 
for the distribution of grain sizes f H3.1|l and the profile 
of optical depth (^21 in the AU Mic disk. Included 
in our analysis are estimates of the sizes and total mass 
of the largest parent bodies undergoing a collisional cas- 
cade. According to our theory, all of the current optical 
to near-IR observations probe grains whose population 
is maintained in steady state and that are still bound — 
most only barely — to the central star. Contributions to 
scattered light from unbound and unequilibrated popula- 
tions of grains are assessed in H3.3l and H3.4I respectively. 

3.1. Equilibrium Size Distribution 

We posit an annulus of parent bodies extending from 
r = tbr — Ar/2 to r — tbr + Ar/2 — the "birth ring" — 
where grains are born from the collisional attrition of 
larger parent bodies. These grains travel on orbits whose 
eccentricities and semi-major axes are continuously mod- 
ified by CPR drag. Grains are removed from the birth 
ring either by CPR drag or by collisions with other 
grains. The goal of this section is to determine the equi- 
librium size distribution dN/ds as a function of r, where 
dN is the vertical column density of grains having sizes 
between s and s + ds. 

The size distribution of dust in debris disks is com- 
monly assumed to be proportional to s~ Qcc = s~ 7 / 2 . This 
scaling is appropriate for grains on low-eccentricity orbits 
whose collisional strengths are independent of size, whose 
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spatial distribution is homogeneous, and whose numbers 
are maintained in a purely collisional equilibrium, as first 
derived by Dohnanyi (1969; see also O'Brien & Green- 
berg 2003; Pan & Sari 2005). In such equilibrium cas- 
cades, as much mass is ground into every size bin as 
is ground out. One presumption behind the cascade is 
that collision times are short enough that the system has 
relaxed into collisional equilibrium. However, Figure ^ 
shows that tcPR and t co i are of the same order for some 
grain sizes. Visible dust may be removed too quickly by 
CPR drag to participate in a purely collisional, equilib- 
rium cascade. 

Instead of making the usual assumption that the size 
distribution is proportional to s~ 7 / 2 for all s, we con- 
struct the following model. We define Sbreak as the radius 
of the grain for which 



icol(sbroak) = ^CPR^break) 



(17) 



at r — tbr. We expect that grains of a given size can 
participate in a standard collisional cascade if they have 
had enough time to collide destructively about once (for 
numerical estimates of the time required for a cascade to 
equilibrate see, e.g., Campo Bagatin et al. 1994 and refer- 
ences therein). For s > s^eak, we expect icPR > t co \ and 
a Dohnanyi-likc size distribution (modified appropriately 
for grains on highly eccentric orbits, i.e., for spatial inho- 
mogeneity). For s < Sbreak, we will se e tha t tppR . < ^col- 
The two regimes are treated in 33.1.21 and 33.1.11 respec- 
tively. We discuss which grain sizes carry the bulk of the 
optical depth and how collision times vary with grain size 
in 33.1.31 These considerations are applied to computing 
Sbreak in 33.1.41 The sizes of the largest parent bodies 
are estimated in 33.1.51 

Readers interested only in our results for the grain size 
distribution can examine equations l|18|) , (|21|l , (|23|l , (|25(1 , 
(j2fj|) . (|2Tjl . l|37jl . and (23), and can skip ahead to 31H1 
where grain dynamics are analyzed. 

3.1.1. Equilibrium Size Distribution for s < Sbreak 
First we define 



dAf 
ds 



dN 
ds 



■ 2nr dr 



(18) 



as the size distribution of grains integrated over the entire 
disk. 

Fresh debris having s < Sbreak continually sprays from 
colliding bodies having s > s^eak- We assume that the 
initial or "injection" spectrum of fresh debris follows a 
power law; the rate at which grains having sizes between 
s and s + ds are injected into the entire disk obeys 



dN 
ds 



(19) 



where C is a constant and the subscript "I" denotes "in- 
jection." Theoretical considerations of mass c onservation 
suggest inject = 3-4 (Greenberg et alJll978(l . while im- 
pact experiments using c entimeter-sized targe ts suggest 
values of ginjcct ~ 3.5-4 ijFuiiwara et al.lll~989i see their 
Figure 3). 

To solve for the steady-state size distribution, we 
equate the injection rate dM /ds\\ to the removal rate 



dAf/ds\n. We show in 33. 1.31 that removal is dominated 
by CPR drag onto the central star: icPRAcoi < 1 for 
s < Sbreak- Then grains having s < Sbreak drag inward 
from the birth ring largely unimpeded by collisions, and 



dN 
ds 



1 



dAf 



*cpr(s) ds 



(20) 



Equating the injection and removal rates yields 



— Cs q ' n]OCt t CPR (s,r pe 

ds 



?*Br) for S < Sbreak ■ 
. . 

The column density dN / ds local to the birth ring is 
proportional to dAf /ds times the fraction of time spent 
inside the birth ring: 



dN 
ds 



x 



dAf Aij-eside | 



BR 



ds i C PR 

s -9i„jcct 



csidc 



CPR 



(22) 



for S < Sbreak , (23) 



where A^esidelcpR is the total time a grain spends in- 
side the birth ring over its CPR-limited lifetime. Figure 
[3 plots At res ide|cPR as a function of s for the parame- 
ters M*/M Q = 1, Ar/r B R = 0.1, and r pcri ,o = t"br- We 
compute this quantity by explicitly tracking the position 
of a grain on a decaying orbit. Evidently, A£ res i(je | cpr 
scales approximately linearly with s for all s. We can un- 
derstand this linear scaling by examining two extremes. 
For s ^ Sbiow, orbits are nearly circular always, and 
At res ide|cPR is merely the time for the grain's orbital ra- 
dius to decay by Ar/2. This time is proportional to s 
since E 1 (see equation Q21)- For s — Sbi ow "C Sbi ow , 
eccentricities eo are large. Over most of a grain's life- 
time, the grain's periastron lies inside the birth ring while 
its apastron lies well outside. Evaluated over intervals 
shorter than £cpr, the fraction of time the grain spends 
inside the birth ring is given by /(e) with r pcr ; set equal 
to Tbr- Applying (|TT)l and i(T5|) . we have 



Ai r , 



'idclcPR' 



oc s 



(X s , 



Ar 

7"BR 



1/2 



l + e 
1 + e 



e_\ 3/5 de 
eo/ e 
(24) 



where the integral in the second row is nearly constant 
with s. 

3.1.2. Equilibrium Size Distribution for s > Sbreak 

For s > Sbreak, a collisional cascade is established be- 
fore CPR drag has time to remove grains. Collisions 
occur primarily in the birth ring since the vertical opti- 
cal depth is greatest there. In the birth ring, the usual 
collisional equilibrium implies 



dN 
ds 



= s 7/2 for s > 



Sbreak • 



BR 



By the same logic that led to (|23|) . 



dAf dN 
ds ds 



f-col 



BR 



At 



for S > Sbreak , 



reside 



(25) 



(26) 
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the regime s < Sbreak? From (|23|) . 
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Fig. 2. — The length of time At rcs id |cpR a grain spends within 
the birth ring (at radii between reR — Ar/2 and rp,R + Ar/2) as 
a function of grain size, for grains whose lifetimes are limited by 
CPR drag. We take r pcri = r B R = 43 AU, Ar/r BR = 0.1, and 
M* = 1M . 



where we have assumed (and show in H3.1.3[) that t co \ < 
tcpR is the appropriate lifetime for grains having s > 
Sbroak- By analogy to Ai resi de|cPR, Ai rosido | co i is the total 
time a grain spends within the birth ring over a collision- 
limited lifetime. 
In the special case of large initial eccentricity eo , 



dN dN 
ds ds 
dN 
ds 



1 



BR /( e o) 

(l-e )- 3/2 . 

BR 



> Sbroak, e f» 1 , (27) 



where we have used the fact that grains having large eo 
and whose lifetimes are limited by i co i < icPR have their 
periastra within the birth ring and their eccentricities 
close to their birth values for nearly all their lives. 

3.1.3. Optical Depth and Collision Times in the Birth 

Ring 

We estimate the sizes of grains that carry the lion's 
share of the optical depth in the birth ring. For 
s > Sbreak, the column density in the birth ring obeys 
dN/ds oc s~ 7 / 2 . The vertical optical depth contributed 
by such grains equals 



/ \ 3^ -1/2 

T_1_(S > Sbreak) ~ S — OC S ' \ 

ds 



(28) 



therefore among grains of size s > Sbreak, the optical 
depth is dominated by grains of size Sbreak- What about 



r±(s < Sbreak) ~ s 3 ^ a s 3 AiresidelcpR ■ ( 29 ) 

For (finject ~ 3.5-4, this quantity either grows or is ap- 
proximately constant with s. We conclude that grains 
of size s ~ Sbreak dominate the total geometric optical 
depth in the birth ring: 

Tj_(sbreak) ~ Tj_ iB R , (30) 

which in combination with J2J , © an d (|10|> implies that 

^col(sbreak) ~ ^c,Br(1 — e 0: brcak) _3 ^ 2 , (31) 
Where e , break = e (s b rcak)- 

We exploit ©, JXOJ , and lf2S J) -lpn jl to estimate t col for 
arbitrary s. For s > Sbreak) 



(32) 



/ s \V2 

t co \{s > Sbreak) ~ ^c,BR ) (1 - e )~ 3/2 , 

\ Sbreak / 

where we have set e = e since most of the grain's lifetime 
is spent with that eccentricity (see M2.4|) . Since in the 
large-s limit t co \ oc s 1 / 2 while icPR oc s, the assumption 
made in H3.1.2I that grains are removed principally by 
collisions for s > Sbreak is asymptotically valid. In Figure 
0] we plot l|3^|) as a fiducial size-dependent collision time, 
replacing Sbreak by Sbi ow to render the curve independent 
of stellar mass- loss rate. This replacement, performed 
solely for Figure ^ accrues only a slight error since we 
find in ij3.1.4l that Sbreak and Sbiow are nearly the same. 

Next we estimate t co \(s < Sbreak) and show that 
tcpn/tcoi < 1 for s < Sbreak, as was assumed in H3.1.1I 
Equation l|30() implies that grains having s < Sbreak are 
predominantly destroyed by Sbreak-sized grains. Then 7 



£col (s < Sbreak) ' 



r2(s)rj_(sbreak) 



' ^col (Sbreak) 



1 — ep.br. 

1 - e 



oak 



3/2 



(33) 



For convenience, we construct the approximate fitting 
formula for t C pR from O , > an( l 03 : 



*CPR(s < Sbreak) ~ ^col(Sbreak) 



Sbreak 



1 - e 



break 



l - eo 



1/2 



(34) 

Combining (|33|l and H34fl . we find that the ratio between 
CPR and collision lifetimes is 



^CPR 



s<s bll 



Sbreak 



l-e 



1 — Co, break 



< i ; (35) 



7 Equation 1331 overestimates t co i because it neglects the fact 
that grains on highly eccentric orbits intercept an optical depth 
parallel to the disk in addition to Tj_ br- This neglect does not 
change our derived scaling relations, but it will change certain nor- 
malizations, e.g., the threshold M* dividing CPR-dominated from 
collision-dominated behavior. We are indebted to Y. Wu for point- 
ing out this omission, which will need to be corrected in future 
work. 
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indeed this ratio vanishes as s approaches Sbiow We have 
established that inequality (|35|1 holds at r = tbr, but in 
fact it holds for all r — r pcrl < tbr, since icPR <x r 2 
while t co i oc r 3 / 2 . Therefore our assumption that bound 
grains having s < Sbreak are removed principally by CPR 
drag is valid. 

3.1.4. Calculating Sbreak 
Equate the collision time (|32[1 to the CPR drag time 



(1 — eo, break) 

to find 



Sbreak 



-3/2 47rc pr BR . 

ic,BR ~ OT p Sbroak-Eieo, break) 



(36) 



1.002s b iow 


= 0.2 /j,m 


for 


M*/M @ 


= 1 


l.OlSblow 


= 0.2 fim 


for 


M*/M© 


= 10 


1-lSblow 


= 0.3 /im 


for 


M*/M 


= 10 2 


2.3sbiow 


= 2 /im 


for 


M*/Mq 


= 10 3 



(37) 

3.1.5. Larger Parent Bodies 

Bodies having s > Sbreak participate in a collisional 
cascade in which the Dohnanyi-like spectrum extends 
from Sbreak up to a maximum size s max . By definition, 
Smax characterizes those grains whose collisional lifetimes 
equal the age of the AU Mic disk; by (|32() . this size sat- 
isfies 



^col(Smax) ^ ^c. 



BR 



Sbrcak 



1/2 



(38) 



where we have dropped the eccentricity-dependent factor 
since grains having s = s ma x S> Sbiow travel on essentially 
circular orbits. Then 



Sbreak 



'c.BR 



10 cm 


for M*/Mq 


= 1 


10 cm 


for M*/M 


= 10 


20 cm 


for M,/M & 


= 10 2 


100 cm 


for M*/M 


= 10 3 






(39) 



Our model has no information on parent bodies hav- 
ing s > s max . While such bodies likely exist, we do not 
know whether they are currently in a constructive (planet 
building) or destructive (debris generating) phase of their 
evolution. It is not justified to extend the size distribu- 
tion to s > s max (e.g., to the kilometer size range) as is 
sometimes done. 

3.2. Physical Implications of Optical Depth Profiles 

The birth ring divides the inner disk from the outer 
disk. In ^3.2.11 we discuss the vertical optical depth in 
the inner disk. In ^3.2.21 we estimate the total disk mass. 
In ^3.2.31 we derive analytically how the vertical optical 
depth scales with radius in the outer disk. 

3.2.1. Inner Disk (r < tbr).' Competition Between 
Collisions and CPR-Driven Accretion 

We showed in H3.1.3I that grains of size s ~ Sbreak 
make the largest contribution to the total optical depth 



at r « tbr. Bound grains having s < Sbreak tend to be 
transported inwards by CPR drag, unimpeded by inter- 
particle collisions. Larger grains tend to be collisionally 
destroyed before reaching the star. How does the vertical 
optical depth in the inner disk compare with the optical 
depth in the birth ring? 

We define CPR-dominated type A disks to be systems 
for which Sbreak — Sbiow Sbiow In such disks, grains 
for which Sbiow < s < Sbreak are numerous, contain a 
significant fraction of the total optical depth in the birth 
ring, and tend to accrete onto the central star before un- 
dergoing a collision. From continuity, the optical depth 
T±(r < tbr) scales approximately as r°: the inner disk 
is "filled in." 

By contrast, in collision-dominated type B disks, 
Sbreak — Sbiow <C s b i ow • The region inside the birth ring is 
virtually empty. The reasons for this are twofold. First, 
the range of sizes of grains that drag in without being 
collisionally destroyed (sbiow < s < Sbreak) is narrow; 
there are not many such grains. Second, because Sbreak 
is so close to Sbiow, Sbreak-sized grains have large initial 
eccentricity. They spend most of their lifetimes having 
fperi.o and only a small portion at r < tbr (see 



As judged from J57)l. if M*/M > 10 2 , type A condi- 
tions would hold for the AU Mic disk and the inner disk 
would be filled in. If M*/M < 10 2 , then type B condi- 
tions would hold and the inner disk would be empty. In 
21 we not only check these assertions by detailed Monte 
Carlo simulations of the AU Mic disk, but also decide 
which case is favored by the observations. 

3.2.2. Total Mass of the Disk 

By (|25|l . most of the mass of the disk is contributed 
by the largest grains (s = s max ), since s 4 dN/ds oc s 1 / 2 . 
We scale from the column density of Sbreak-sized grains 
in the birth ring, 



Tj_,BR ~ 7TS brcak 



dN 
ds 



to estimate the column density of s E 
the birth ring, 



(40) 

-sized grains in 



dN 
ds 



,BR 



7rS break 



Sbreak 



-7/2 



(41) 



The total number of such grains is their column density 
multiplied by the area of the birth ring 27rrBRAr, since 
Smax-sized grains undergo a collision long before dragging 
inwards. Multiplying this total number by the mass of a 
single grain yields the mass of the disk, 



M„ 



8tt 



-p»"BRTJ-,BRSbrcak 



Sbreak 



1/2 



Ar 



(42) 



'O.OlMff 



Ar/rgR \ / ti,br 
) \ 4x 10- 3 



"break 



0.1 / V4 x 10- 3 / \0.2pim / 

where we have used and 139f) . and have normalized 
Ar/rBR and ti.br to values that yield good fits to ob- 
servations, as described in The steady comminu- 
tion of ~O.O1M ~ 1 lunar mass's worth of decimeter- 
sized bodies into micron-sized particles does not seem 



an unduly heavy burden for the AU Mic system to 
bear. The solar system is thought to have somehow shed 
~10Mq of rock and ice near ~30 AU over an uncertain 
timescale of 10-1000 Myr during its "clean-up" phase 
^Goldreich. Lithwir.k. fc Sarill20f)l . 

3.2.3. Outer Disk (r > tbr).' Barely Bound Grains 
Grains created at r = tbr = 43 AU and having 



Sblo 



< 



Sblo 



having periastron distances 



pen 



?"br (see ^2.3|l . Here 



we show that such barely bound grains establish an outer 
disk at r > tbr whose vertical optical depth scales ap- 
proximately as r~ 5 / 2 for CPR-dominated type A disks, 
and as r -3 / 2 for collision-dominated type B disks. The 
contribution of unbound grains having s < Sb iow r elative 
to that of barely bound grains is assessed in H3.3I 

The outer disk comprises grains having sizes slightly 
greater than Sbiow since only those grains have substan- 
tial eccentricities (see equations 0). We refer to such 
barely bound grains as having sizes s — s D iow.+ - The op- 
tical depth in the outer disk should scale approximately 
as their column density iVbi ow .+ : 



T|_ OC iVblow,+ °C -- 



1 dA/i 



blow,+ 



, * (43) 
r dr 

where A/blow, + is the total number of such grains in the 
entire disk (sec |18|). By the chain rule, 



blow, 



blow,+ 



de 



dr 



de dr 



(44) 



We determine de/ dr by making the approximation that 
at any instant, all grains are located at their apastra: 



(45) 
(46) 



' apo l ' pen 

r - r B R 



This approximation should be good for the highly eccen- 
tric orbits of the outer disk. Differentiating H46[) . we find 
for r ^> rBR that 



2r B R 



de 



dr (r + trr) 2 



1 



(47) 



The remaining factor in 144|) . dA/biow.+/^e, is deter- 
mined by the size distribution of barely bound grains. 
This differs between type A and type B disks, as seen 
below. The size distribution determines the initial shape 
of the eccentricity distribution. The eccentricity distri- 
bution is altered over time by CPR drag. Consider grains 
having identical initial eccentricities eo created at a con- 
stant rate Af. In steady state, CPR drag transports a 
constant number of particles per time through eccentric- 
ity space, (dAf /de){de/dt) ~ Af. Hence, 



dAf 



— — oc . , 

de \ at 



de 



,3/5 



x 



(1 - e 2 ) 3 / 2 



for e < eo ■ 



(48) 



where we have used ljT5|) . Equation (|4*8|) implies that 
nearly all barely bound grains have e « eo . Therefore the 



CPR-evolved eccentricity distribution closely resembles 
the initial eccentricity distribution: 



dAf dAf 
de deo 



(49) 



occupy initially highly eccentric orbits dAf/d, 



Now we address the size distribution that determines 
dAf /deo. For type A disks, for which Sbrcak — Sbiow 
Sbiow j the size distribution for barely bound grains obeys 



s cx s 9injcct icpR(s) (see |2T]). For type B disks, 
Sbroak- Sbiow < Sbiow! outer-disk grains for which s b i ow < 
s < Sbrcak ar e outnumbered by grains having s > Sbreak- 
The relevant size distribution in the outer regions of type 
B disks is therefore the one appropriate for s > Sbreak^ 
dAf/ds oc s~ 7 / 2 (l-e )~ 3/2 (see E3). 8 We evaluate these 
distributions for s = Sbiow,+ : 



dAf hlow ,+ ^ £ (eo) k (1 _ eo) _ 1/2 for Type A digks ( 



ds 



(50a) 



dA/biow,+ _ eo) _ 3/2 forTypeBdiskSi (50b) 
ds 

using equations $T2$ and l(TB|l . Now dAf/ de ~ dAf/ deo = 
(dAf /ds)(ds/deo). Since equations © imply that ds/deo 
is approximately constant for eo ~ 1 (i.e., for s = 
Sbiow, +), and since r » r apo » 2r B R/(l - e ), 



dA/ " blow ^ + -x (1 - e )- 1/2 oc r 1 ' 2 for Type A disks, 



de 



(51a) 



blow ' + oc (1 - e )- 3/2 cx r 3 / 2 for Type B disks . 



de 



(51b) 



Combine |@3J|, igU, 1071), and JHH to obtain 

ti oc r _5/2 for Type A disks , (52a) 



tj_ oc r 



-3/2 



for Type B disks . 



(52b) 



Note that these scaling relations cannot be obtained by 
merely assuming that the disk-integrated size distribu- 
tion obeys the usual Dohnanyi relation dAf /ds oc s -7 / 2 . 
For type B disks, for example, the key modification arises 
from the prolonging of the collisional lifetime due to SWR 
pressure (i.e., the factor of (1 — eo)~ 3 / 2 in equation |51bj ). 
Moreover, the scalings are robust against uncertainties in 
the size distribution; they do not depend explicitly on ei- 
ther q ce or (finject- We verify these scalings by numerical 
experiments in ^ 

3.3. Unbound Grains ((3 > 1/2) 

The rapid expulsion of unbound (/? > 1/2) grains com- 
pared to the longer, CPR-driven orbital decay of barely 
bound (1/2 — /3 <C 1) grains suggests that in steady state, 

8 This statement for type B disks is only valid not too far from 
the birth ring. As r — > oo, the only bound grains that are present 
have sizes between Sbiow and Sbrcak- These obey tj_ oc r~ 5 / 2 , 
just as they do for type A disks. Therefore for type B disks, the 
t± oc r~ 3 / 2 scaling derived in the main text eventually gives way 



to tj_ cx r 



-5/2 
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unbound grains contribute little to the surface brightness 
of the outer disk. On the other hand, unbound grain 
velocities are nearly constant with radius — for j3 w 1, 
velocities are approximately equal to their (circular, Ke- 
plerian) birth velocities. As a result, the optical depth 
of unbound grains should roughly obey rj_ iU b oc r _1 and 
should eventually exceed, at some "cross-over radius," 
the optical depth of barely bound grains, which scales 



as ri.bb oc r 



-5/2 



in type A disks and as rj^bb oc r 



-3/2 



in type B disks f ^3.2.3|) . In this section, we estimate 
the value of the cross-over radius, r cross , and show in the 
case of AU Mic that it lies outside the scope of current 
observations. 

Consider a type A disk. In the birth ring, the optical 
depth of barely bound grains exceeds that of unbound 
grains by 



T"L,bb/Tl,ub 



dN 

ds 



s 2 ds 



Sblow dN 

min(sv) ^ 



s 2 ds . 



( 53 ) 

The smallest unbound grain of interest is the smallest 
grain for which Q sca t ~ 1: min(sy) « 0.1 /im. Since 
■?\L,bb oc r~ 5 / 2 while TL iU b oc r — 1 , the cross-over radius is 



^ 2/3 
Tl,bb \ 

Across ~ 1 ) r B R 

TL.ub 



(54) 



By m, 



dN 
ds 



~ Ds 9ln j ctt At res i(je I CPR, (55) 

Sblow<s<2s b iow 

in the birth ring, where D is a constant. By the same 
logic that led to (J23J, 



dN 
ds 

S<Sblo 

in the birth ring, where 



D*-*»J"*tbiow(fl) 



(56) 



iblow(s) 



Vr B R.Ar 



90 



Sblo 



1/2 



Ar 

— yr (57) 
r-BR 



is the time for an unbound grain to leave the birth ring, 
and f^BR = yjGM*/r^. Numerical evaluation of the 
integrals in l|53|l reveals that 



900 AU 
200 AU 



(58) 



for M*/M Q = 10 2 
for NU/Mq = 10 3 . 

If instead type B conditions apply for the AU Mic 
disk, then by considerations analogous to those above, 
Across 3> 10 3 AU. We compare r cross with the maximum 
radius probed by current observations — approximately 
200 AU — to conclude that under type A conditions, un- 
bound grains contribute at most marginally to the cur- 
rently observed surface brightness. Under type B condi- 
tions, they contribute negligibly. 

3.4. Unequilibrated Grains (Sblow < s < s age j 

Grains on extremely eccentric orbits may have colli- 
sional and CPR lifetimes that exceed the age of the sys- 
tem. Can such grains, whose numbers cannot be assessed 



within our steady-state model, contribute significantly 
to the observed surface brightness of the AU Mic disk? 
Such "unequilibrated grains" have sizes between Sbiow 
and Saec, where 



iago — niin[icPR(s a go), ^col(Sage)] 

defines s ag0 . We apply i|12[l and (|33|l to find that 



(59) 



Sblo 



1.006 
1.002 
1.00003 
1.000002 



(60) 



for M*/M Q = 1 

for M*/M@ = 10 

for M*/M@ = 10 2 

for M*/M s = 10 3 . 

Unequilibrated grains occupy such a narrow range of 
sizes that they seem unlikely to contribute much to the 
total optical depth. We ignore the unequilibrated popu- 
lation for the remainder of this paper. 

4. MONTE CARLO MODELLING 
4.1. Procedure 

To test several of the analytic results derived in <JU we 
model the AU Mic disk by means of a Monte Carlo sim- 
ulation. We calculate the geometric optical depth T±(r), 
edge-on surface brightness profile SB(6), and spectral flux 
density F u and then we compare to observations. The in- 
put parameters are the stellar mass- loss rate M» , optical 
depth in the center of the birth ring t_l,br = Tj_(?"br)j 
and width of the birth ring Ar. 

We lay down a number J of dust particles around the 
central star in a two-dimensional plane. Each particle's 
radial coordinate r and azimuth ip are determined by the 
particle's semi-major axis a, eccentricity e, true anomaly 
4>t, and longitude of periastron uj. Since our model disk 
is axisymmetric, lu is drawn as a uniform deviate from 
to 2tt. 

The birth distributions of the remaining orbital ele- 
ments depend on the distribution of grain sizes s. Or- 
bital elements subsequently evolve from their birth values 
by CPR drag. The degree of evolution depends on the 
age of the grain (t) relative to the CPR lifetime (£cpr! 
equation |12j ) and collision lifetime (t co \; equations |32j 
and |33j). By definition, Sbrcak is the grain size for which 
^cpr — ^coi, and s age is the smallest grain size for which 
min(t C PR, icoi) = £ ag c- Only grains having s > s agc can 
be removed over the age of the system. For s > Sbreak, 
collisions are more important than CPR drag in remov- 
ing grains (t co \ < £cpr)j so we draw the ages t of such 
particles as a uniform deviate from to t co \(s). Grains of 
size s < Sbrcak are removed by CPR drag (tcPR < *col), 
so we draw t for these grains as a uniform deviate from 
to £ C pr(s). 

The evolved eccentricity e depends on the initial eccen- 
tricity eo and the age of the particle, implicitly according 
to 



t 



= 1 - 



E(e) fl + 



eo 



-) 

eo/ 



8/5 



(61) 



*cpr E(e ) \ 1 

where E(e) is defined by (|14|) . The evolved semi-major 
axis is given by 



e 

ao 1 — 



4/5 



1 



1-e 2 



(62) 
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l)Wvatt fc Whipplel Il950j) . where the initial semi-major 
axis clq is derived from the initial periastron r por i,o, which 
we draw in the following way. In our model, all grains are 
born in the birth ring, an annulus of width Ar centered 
at tbr- At birth, a grain is located at the periastron of 
its osculating orbit — an orbit rendered eccentric by SWR 
forces ( M2.3[l . We draw r per i,o from a uniform distribution 
of width Ar centered at tbk- 

As mentioned, the birth distributions of eccentricities 
and semi-major axes depend on the distribution of grain 
sizes s. Although particles in our simulation are born 
only in the birth ring, their steady-state population may 
occupy all space, so we must draw s from the global, 
disk- integrated size distribution dM /ds. We apply re- 
sults from H3 . II made more precise for our Monte Carlo 
calculation. From l|21|l . 



dM 



CX S 9inicct t C PR(s,rpcri = Tbr) for S agc < S < Sbreak , 

(63) 



ds 

and from i|26|) . 



dM 



t, 



col 



for S > Sbreak • 



(64) 



ds A£ rcs ide| co i 

We approximate t co \ using l|32|) . and we evaluate 
Ai res ide|coi numerically. We take the distribution de- 
scribed by (|fJ3*|) and l|6^|l to be non-zero only for s > s agc 9 
and to be continuous across Sbreak; furthermore, we trun- 
cate the distribution at s max = 500sbi ow because of com- 
putational limitations. We fix <7inj C ct = 4. The distribu- 
tion dM I ds is plotted in Figure |3] For each s drawn, eo 
and ao are calculated using equations JJjJ. For J = 10 9 , 
one Monte Carlo simulation takes 8 hr to complete on a 
1.33 GHz PowerPC G4 processor. 

That the distribution of mean anomalies <fi m is uniform 
determines the distribution of true anoma lies <pt via Ke- 
pler's equation l|Murrav fc Dermottlll999lk 



<Pm = <Pc - e sin e , 
where <p is the eccentric anomaly: 

1/2 



. <f>t A + e 

tan — = 

2 V 1-e 



tan ■ 



(65) 



(66) 



Knowing a, e and <pt for each particle determines its ra- 
dial distance from the star: 



a(l-e 2 ) 



(67) 



1 + e cos (f>t 

4.2. Products of the Monte Carlo Calculation 

Having laid down J particles of various sizes, we output 
the following: 

1. The geometric vertical optical depth Tj_(r). We 
first calculate this quantity in relative units by 
summing the geometric cross sections of particles 
in a given annulus, and dividing the resultant sum 
by the area of that annulus. We then normalize this 
result by matching the input t^br to the model's 
relative optical depth at tbr. 

9 For M,/Mq = 1, Sago > screak so on ly IG4I is relevant. 
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Fig. 3. — Disk-integrated grain size distributions dj\f/ds. Dashed 
and triple-dot-dashed curves correspond respectively to M, / Mq = 
{10, 10 3 }. A dotted line proportional to the Dohnanyi scaling 
s — 9oc = s — 7/2 ; s overplotted for reference. As s approaches Sblow> 
the population of grains rises significantly above what one would 
expect from a pure Dohnanyi size spectrum. The deviations from a 
Dohnanyi spectrum, which differ under type A (e.g., M„ = 10 3 Afo) 
and type B (e.g., M* = IOMq) conditions, arc critic al for under- 
standing how tj_ scales with r in the outer disk; see H3.2.3I 

2. The surface brightness of the disk observed edge- 
on at y-band (Hubble Space Telescope's F606W) 
and iJ-band wavelengths, as a function of projected 
stellar separation b: 



SB(6, A) 



^% Qscat (A, s) P(8, A, s) its 2 ^{r, s) di ds , 
47rr z ds 



(68) 



where £ = ±vV 2 — b 2 measures distance along our 
line of sight. The stellar flux incident on a grain 
in the wave band of interest is AL> j »/47rr 2 , the 
cross section for scattering is Qscat tts 2 , the scat- 
tering angle between the star, grain, and observer 
is 9 = tan _1 (£/6), and the relative power scattered 
per steradian is P (normalized so that its integral 
over all solid angle equals unity). We use Mie the- 
ory to calculate Q scat and P, adopting the optical 
constants of pure water ice ijWarrenl 11984'). The 
volumetric number density of grains n is found by 
dividing the geometric vertical optical depth by the 
height of the disk: 



2 dn 
ns -(r,s) 



1 dr± 
h(r) ds 



(r, s) . 



(69) 



The radial height profile h(r) is derived from the 
projected disk height h(b), which roughly follows a 
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broken power law that changes slope around b 
r BR = 43 AU: 



Kb) 



'BR 



(V^br)" 1 
(V^br)" 2 



for b < r B R 
for b > r B R . 



(70) 



We set h B R = 2.75 AU. iKrist et all (|2005|) fit sepa- 
rate broken power laws to the northwest and south- 
east extensions of the disk and obtain r\\ « and 
r]2 ~ 1-2; see their Figure 7. For simplicity, we 
take 771 = and 772 = 1. These values characterize 
an inner disk that is empty and seen in projection 
(or that has constant height), and an outer disk in 
which grains have constant inclination dispersion. 
We adopt a radial profile h{r) identical to h(b) as 
given in l|7UI) with b replaced by r. We divide the 
F606W profile by the H profile to obtain a V — H 
color profile. 

3. The spectral flux density F v : 



F v = —B V (T*) nRl 

+ Jj ' B v (T(r, s))Q cmis (X,s)7rs 2 "J^O", s) 2irrdr ds , 

(71) 

where B v (T) is the Planck function. We model the 
emissivity of the dust as a broken power law: 



Qe 



1 

2-ks/X 



for 2irs > X 
otherwise , 



(72) 



in approximate ag reement with the model of 
iChiang et ail $Mr$) for ice-mantled silicate grains. 
Since the peak wavelength of emission from AU Mic 
is about 1 /xm, all bound grains are large enough to 
absorb most of the incident stellar flux. We solve 



ITS = 4-7T 



(X,s)irs 2 dX (73) 



for the temperature T specific to a given grain size. 
4.3. Results 

By experimenting with various values of tj^br, Ar, 
and M* <G (1, 10, 10 2 , 10 3 )M Q , we find that r ±3R = 
4x 10~ 3 , Ar/r BR = 0.1, and M* £ (1, 1O)M yield theo- 
retical surface brightness profiles and spectra that agree 
encouragingly well with observations. Our preferred val- 
ues for tj^br and Ar are likely uniquely determined to 
within factors of a few; SB(6 = 7'br) oc tj^brV Ar while 
F v cx Ti,BRAr. The close resemblance of the models 
for which M* < 10M Q means that we cannot do better 
than place an upper limit on M^/Mq of ~10. In what 
follows, we present results for our preferred input param- 
eters, in addition to models for which M* £ (10 2 , 10 3 )AfQ 
to study the effect of varying M* alone. 

Figure 0] displays geometric optical depth profiles 
Ti(r). As expected from our analysis in H3.2I disks 




10 100 1000 
r (AU) 



Fig. 4. — Vertical optical depth profiles computed from 
our Monte Carlo simulations. Solid, dashed, dot-dashed, and 
triple-dot— dashed curves correspond respectively to M, /M© = 
{1, 10, 10 2 , 10 3 }. Values for t x ,br = 0.004 and Ar/r BR = 0.1 are 
held fixed for all models. The two types of disks, CPR-dominatcd 
type A disks and collision-dominated type B disks, may be dis- 
tinguished. The inner disk for M» = IMq is completely empty 
because s broak < s ago . 



separate into two types, CPR-dominated type A and 
collision-dominated type B. As M* increases, disk be- 
havior changes from type B to Type A. One consequence 
is that the inner disk becomes increasingly filled in. Fur- 
thermore, for type A disks, we expect tj_ oc r~ 5 / 2 at 
r 3> rBR; this behavior is indeed evident for M^/Mq £ 
(10 2 ,10 3 ). For type B disks, we expect t± oc r~ 3 / 2 at 
r ^S> rBR; the models for which M»/M Q £ (1, 10) exhibit 
this scaling. 

We compare our theoretical surface brightness profiles 
SB(&) with data from K05 in Figures El and El for the 
cases of high M* and low M* , respectively. For all disk 
models, there is a significant contribution to the surface 
brightness at b < r B R. from starlight that is forward scat- 
tered by grains located within the half of the birth ring 
nearest the observer. As M* increases, the inner disk 
fills in and the surface brightness at b < rBR increases. 
In comparison, the surface brightness at b 3> rBR de- 
creases with increasing M* , reflecting the transition from 
the tj_ oc r~ 3 / 2 scaling of type B disks to the tl oc r~ 5 / 2 
scaling of type A disks. Examination of either the in- 
ner or outer disk profiles reveals that models for which 
M*/Mq £ (1, 10) fit the data better than do models for 
which is higher. Discrepancies between these low- 
M* models and the observations are less than a factor of 
2. They might arise in part from our use of a scatter- 
ing phase function (P) appropriate for idealized spherical 
grains of pure water ice. 

In Figure we plot V — H colors. For low values of 
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10 100 

Projected Separation b (AU) 

Fig. 5. — Theoretical type A and observed surface brightness pro- 
files. Thin lines are data (for northwest and southeast extensions of 
the disk) from K05. dot-dashed and triple-dot-dashed curves cor- 
respond respectively to M*/Mq = {10 2 , 10 3 }. The vertical dotted 
line corresponds to = 43 AU, the radius of the birth ring con- 
taining dust-producing parent bodies. The inset scaling of ft -9 / 2 
is derived from the rule of thumb that at large b, SB oc 6 7— 17-1 for 
h oc r* 1 and tj_ oc r 7 . According to our theory for type A disks, 
r\ = \ and 7 = —5/2. 




10 100 

Projected Separation b (AU) 

Fig. 6. — Same as Figure 151 but for type B disk models. Thick 
solid and dashed curves correspond respectively to M* / Mq = 
{1,10}. The vertical dotted line corresponds to tbr = 43 AU, 
the radius of the birth ring containing dust-producing parent bod- 
ies. The inset scaling of 6~ 7 / 2 is derived from the rule of thumb 
that at large b, SB <x b' y ~ v ~ 1 for h oar 71 and t±_ k r 7 . According 
to our theory for type B disks, r\ = 1 and 7 = —3/2. Stellar mass- 
loss rates of I-IOMq yield surface brightness profiles that agree 
better with the data than those derived from rates of lO 2 -lO 3 A/0; 
contrast with Figure 151 



M*, the inner disk is largely evacuated and so there is 
little variation in color with b for b < t^r- The outer 
disk becomes progressively bluer with b as ever smaller 
(still bound) grains are probed. All of these trends are 
in agreement with observations of disk color (S. Metchev 
2005, private communication; M. Fitzgerald 2005, pri- 
vate communication). 

In Figure we plot our theore tical spec t ra to- 
geth er with flux meas urements from ILiu et all (|2004T) 
and IChen et alJ ((2005) . The filled inner disks of high- 
M*, type A models produce too much emission at mid- 
infrared wavelengths to compare well with observations. 
As was our conclusion from studying the surface bright- 
ness profile in reflected starlight, the thermal emission 
spectra point to stellar mass-loss rates of < 10Af Q . 

5. SUMMARY AND DIRECTIONS FOR FUTURE WORK 

We have constructed a theory to explain the observed 
optical surface brightness profile and infrared emission 
spectrum of the debris disk encircling AU Mic. In 
our theory, the slope of surface brightness versus pro- 
jected radius b changes abruptly at b = 43 AU because 
a birth ring of planetesimals exists at stellocentric ra- 
dius r = tbr = 43 AU. This ring has a full radial width 
Ar ~ OTtbr and a vertical, geometric optical depth of 
Tl.br ~ 0.004. The parent bodies in the ring have sizes 
Smax ~ 10 cm and a total mass of M max ~ 0.01M ffi . Col- 
lisional attrition of parent bodies generates micron-sized 
grains that scatter starlight at optical wavelengths. The 
population of visible grains is maintained in steady state: 
production by colliding parent bodies balances removal 



by grain-grain collisions and removal by corpuscular and 
Poynting-Robertson (CPR) drag. The timescales over 
which removal of visible grains occurs can be orders of 
magnitude shorter than the age of the system (~12 Myr), 
ensuring steady state. 

Collisions between parent bodies initiate a collisional 
cascade that extends downward in particle size by sev- 
eral orders of magnitude. Grains having sizes s < Sbiow ~ 
0.2 jitm are expelled from the system by stellar wind and 
radiation (SWR) pressure and contribute negligibly to 
the observed optical emission. Instead, barely bound 
grains, having sizes just larger than Sbiow and which 
occupy highly eccentric orbits, make the dominant con- 
tribution to the surface brightness in the outer disk at 
r > ^br- The number of such grains rises more steeply 
than would be expected from a pure Dohnanyi size spec- 
trum as s approaches Sbiow from above, because grains on 
high-eccentricity orbits have prolonged lifetimes against 
CPR drag and collisions. The structure of the outer disk 
depends on whether these smallest of bound grains are 
removed principally by CPR drag (type A conditions) or 
by destructive grain-grain collisions (type B conditions). 
As the luminosity and/or mass-loss rate of the central 
star increases, disk behavior grades from type B to type 
A. As the number of parent bodies in the birth ring in- 
creases, collision rates increase and disk behavior changes 
from type A to type B. In the outer reaches of type A 
disks, the vertical optical depth scales approximately as 



ti cx r 



-5/2 



Under type B conditions, t± cx r 



-3/2 



(but 
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Fig. 7. — Color profile (F60GW — H) computed using our 
theoretical Monte Carlo model. Solid, dashed, dot-dashed, and 
triple-dot— dashed curves correspond respectively to M*/Mq = 
{1, 10, 10 2 , 10 3 }. The vertical dotted line corresponds to reR = 
43 AU, the radius of the birth ring containing dust-producing par- 
ent bodies; the horizontal dotted line is the star's color. The outer 
disk is expected to be bluer in scattered near-infrared light than 
the inner disk. 



see footnote [SJ). We have derived these scaling relations 
analytically and have verified them by Monte Carlo sim- 
ulations. 

The inner regions at r < tbr are populated by grains 
that survive long enough before suffering destructive col- 
lisions that their periastron distances diminish apprecia- 
bly by CPR drag. In type A disks, a significant fraction 
of grains born in the birth ring meet this criterion, so the 
inner disk is characterized by the same vertical optical 
depth that characterizes the birth ring. By contrast, un- 
der type B conditions, the inner disk is practically empty. 

In the case of AU Mic, type B conditions prevail. By 
fitting simultaneously both the surface brightness profile 
and the thermal emission spectrum, we not only uniquely 
determine the vertical optical depth and radial width 
of the birth ring (see the values cited above), but also 
constrain the stellar mass-loss rate M* to be < lOM©. 
According to our theory, the inner disk of AU Mic at 
t < tbk is empty. The observed surface brightness at 
b < fBR is not zero because we are observing the disk 
edge-on. The primary contribution to the surface bright- 
ness at b < tbr arises from starlight that is forward scat- 
tered by grains in the birth ring. 

Our theory states that the observed structure of the 
AU Mic disk reflects processes that are balanced in 
steady state. Equilibrium is likely since the timescales 
over which collisions and CPR drag operate, even in 
the rarefied outer disk, are shorter than the age of the 
system. The outer disk does not comprise "primordial" 



Fig. 8. — Spectra computed using our theoretical Monte Carlo 
model. Solid, dashed, dot-dashed, and triple-dot-dashed curves 
correspond resp ectively to M*/Mq = {1, 1 0, 10 2 , 10 3 }. Data from 
ILiu et all 120041) (circles and diamond) and lChen et alj (12003') (tri- 
angles) are overlaid; the open circle is a possible detection, and the 
diamond represents an upper limit. The discrepancy between the- 
ory and observation at A £3 12 fim is due to our use of a blackbody 
spectrum for the central M dwarf; better agreement can be had by 
employing more realistic stellar atmosphere models (e.g., Allard ct 
al. 2001). The discrepancy between the data and the low-M* mod- 
els at the longest wavelengths is due at least in part to the fact that 
the maximum grain size in our simulations is 500sblow 'C s max ; 
were we to increase the number of particles J in our simulation, 
the discrepancy would be reduced. 



grains left behind from a now-evaporated gaseous disk, as 
has been speculated previously. Nor is the manifestation 
of the debris disk phenomenon in AU Mic the outcome 
of a recent cataclysm that has not yet equilibrated. That 
our required parent body mass is modest (equation |42j ) 
supports our contention that the AU Mic disk is in steady 
state. 

As our paper was being completed, we became aware of 
an ind ependent study of AU Mic bv lAueereau fc Beustl 
( 2006). These authors find by inverting the observed sur- 
face brightness profile that the underlying vertical opti- 
cal depth of the AU Mic disk peaks near 35 AU. It is 
heartening that their conclusion is consonant with one 
of ours, derived as they are using complementary ap- 
proaches: detailed data-fitting procedures versus physi- 
cal reasoning to understand dust dynamics under general 
circumstances. 

The ring of parent bodies at tbr = 43 AU that we en- 
vision encircling AU Mic presents a youthful analogue 
to the Solar System's Kuiper belt (see the Protostars 
and Planets V review by Chiang et al. 2006). The spa- 
tial dimensions of these systems are remarkably similar: 
The Classical Kuiper belt, containing those planetesimals 
thought to have formed in situ, extends in heliocentric 
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distance from ~40 AU to ^47 AU (e.g., Trujillo & Brown 
2001). 

We conclude by pointing out directions for future work 
on AU Mic and other debris disks. 

1. Disk thickness. — By assuming that we are viewing 
the AU Mic disk perfectly edge-on, we estimated 
a full disk height of /ibr ~ 2.75 AU. The corre- 
sponding opening angle is /ibrAbr ~ 4°. While 
our model of a type B disk succeeds in reproduc- 
ing the observed scaling behavior of disk height 
(h oc 6° in the inner disk and h oc b 1 in the outer 
disk) , we have not explained what sets the normal- 
ization. Dissipative grain-grain collisions would be 
expected to damp the inclination dispersion and to 
reduce /ibr to values orders of magnitude smaller 
than our inferred value. 

2. Application to other systems. — The debris disk sur- 
rounding the A-st ar 8 Pict oris closely resembles 
the AU Mic disk l|Lml I2004D . The surface bright- 
ness profile abruptly changes slope at b ss 100 AU, 
from SB oc b~ 2A to SB oc 6~ 4 (Kalas & Jewitt 
1995). Moreover, the vertical scale height h scales 
with b the same way that it does in the AU Mic 
disk. 10 Recently, another analogue to the AU Mic 
disk has been discovered: F-star HD 139664 hosts 
a debris disk whose surface brigh tness profile ex - 
hibits a sharp break at ~90 AU ijKalas et alJ l2006) . 
The theory we have laid out for AU Mic might find 
ready application to these other systems. 

3. Uniqueness of AU Mic among M dwarfs. — The 
pioneering Keck survey conducted by PJL05 at 
A = 11.7 /im reveals that AU Mic is distinguished 
among their sample of nine M dwarfs having ages of 
10-500 Myr in emitting an infrared excess. Why? 
Do the other M dwarfs not possess disks? As M 
dwarfs constitute the most numerous stars in the 
universe, understanding why AU Mic might be ex- 
ceptional will help to determine the prevalence of 
planetary systems. Many of the M dwarfs surveyed 
by PJL05 may simply be much older than AU Mic; 
their parent body populations may have suffered 
near complete comminution. 

4. Ubiquity of rings. — That parent bodies are con- 
fined to a ring centered at 43 AU in the AU 
Mic system calls for explanation. Ring morpholo- 
gies are so common — witness the examples of HR 
4796A (Schneider et al. 1999), e Eridani (Greaves 
ct al. 1998), Fomalhaut (Kalas et al. 2005), and 
even the Kuiper Belt — that the "debris disk phe- 
nomenon" might well be more precisely termed the 
"debris ring phenomenon." While regions interior 
to rings might have been purged of material by 
planets, the physical processes that determine the 
outer edges of rings remain unclear. Ideas proposed 
by Takeuchi & Artymowicz (2001) and Klahr & Lin 

10 That the color of Pic's outer disk is red rather than blue 
(by contrast to the case of AU Mic) could be a consequence of the 
particular grain size required for blow-out in the j3 Pic system, since 
for certain grain sizes and compositions, Qscat(A,s) can actually 
increase with increasing wavelength. See Bohren & Huffman (1983) 
for a discussion of this phenomenon of "blueing." 



(2005) for how interactions between solids and gas 
can concentrate planetesimals into rings might be 
relevant. That planetary systems have sharp outer 
edges suggests that planetesimal formation is not a 
continuous function of disk properties; rather, the 
formation of planets may require disk properties to 
meet threshold conditions (e.g., Youdin 2004). 
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TABLE 1 

Stellar Wind and Radiation Parameters in the AU Mic System 51 



M,/Mq 


fsWR 


PCPK 


«blow(A"ll) b 


*CPR(s = l-ls b low)(yr) c 


*cprO 


i = l-5s blow )(yr) c 


<CPr(s = 15s blow )(yr) c 


1 


2.0 


5.0 


0.23 


9.5 x 10 6 




4.9 x 10 6 


1.9 X 10 7 


10 


2.0 


32 


0.23 


1.5 x 10 6 




7.8 x 10 5 


3.1 x 10 6 


10 2 


2.4 


3.0 x 10 2 


0.28 


1.9 x 10 5 




9.9 x 10 4 


3.9 x 10 5 


10 3 


6.4 


3.0 x 10 3 


0.74 


5.1 x 10 4 




2.6 x 10 4 


1.0 x 10 s 


a Assumcs Q ra d — 


2, Qwind = 


1, p = 2g cm -3 


w wind = 450 km s -1 , L, = 


0.1L e , 


M, = 0.5M o , M Q 


= 2 X lCr 14 Af yr -1 , and 



i-pori = r B R = 43 AU. 
^The blow-out radius is such that {3 — 1/2. 

c Calculated for e = e . For s/s blow = {1.1, 1.5, 15}, e = {0.83, 0.50, 0.034} by equations J9) and E(e ) = {7.2, 2.7, 1.1} by equation 
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